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Abstract 

Complete energy spectrum is obtained for the quantum mechanical problem of N 
one dimensional equal mass particles interacting via potential 

N i 

1 a 



V(x 1 ,x 2 ,...,x N )=gJ2 ( x ._ x .y 



Further, it is shown that scattering configuration, characterized by initial momenta 
Pi(i = 1,2,..., N) goes over into a final configuration characterized uniquely by the final 
momenta p[ with p[ = piy + i-i. 
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In recent years, the Calogero - Sutherland (CS) type of N-body problems in one 
dimension have received considerable attention in the literature [1,2,3,4]. It is believed 
that the CS model with inverse square interaction provides an example of an ideal gas in 
one dimension with fractional statistics [5]. Besides, these models are related to (1+1) - 
dimensional conformal field theory, random matrices as well as host of other things [6]. 
Inspired by these successes, it is of considerable interest to discover new exactly solvable 
N-body problems. 

The purpose of this note is to present one such example. In particular I show that 
the N-body problem with equal mass in 1-dimension characterized by (h = 2m = 1, g > 
-1/2, a > 0) 

N Q2 n 
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is exactly solvable. The interesting point about this model is that unlike most other 
exactly solvable models, it has both bound state and scattering solutions. In particular 
I show that the complete bound state spectrum ( in the center-of-mass frame) is given 
by the formula 

where 

, N(N-l) N(N + 1) 3 1 f- — 

For positive energy one has only scattering states. I show that a scattering configura- 
tion, characterized by initial momenta pi(i = 1, 2, N) goes over into final configuration 
characterized uniquely by the final momenta p[ with 

Pi = PN+l-i (4) 
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However, unlike the pure inverse square scattering case (a = 0), in our case the phase 
shift is energy dependent. Thus, as in other integrable cases, in our case too the scattering 
problem reduces to a sequence of 2-body processes. 

Finally, a la Sutherland [3] , I also solve a slightly diferent variant of the Hamiltonian 
(1) with —a/ \jYli<j( x i ~ x j) 2 being replaced by an external potential —a/ a/X^ x ] an d 
obtain exact expressions for the ground state energy eigenvalues and eigenf unctions. 

Consider the Hamiltonian as given by eq. (1). We need to solve the eigenvalue 
equation 

Hip = Eip (5) 

where ip is a translation invariant eigenfunction. Note that our Hamiltonian is very simi- 
lar to the classic Calogero Hamiltonian (see eq. (2.1) of his paper [2]) except that whereas 
he has a pairwise quadratic potntial, we have a "N-body" potential -, ~ a How- 

ever, we shall see that many of the key steps are very similar in the two cases and hence 
we avoid giving most steps which are already contained there [2]. Without any loss of 
generality, a la [2] we also restrict our attention to the sector of the configuration space 
corresponding to a definite ordering of particles, say 

Xi > x i+ i, i = 1,2, N- 1. (6) 

A la [2] it is clear that the normalizable solutions of eq. (5) (with H being given by eq. 
(1)) can be cast in the form 

ij(x) = Z a + 1 / 2 <j ) (r)P k (x) (7) 

where a is defined in eq.(3) while Z and r are given by 

1 N 

Z = U^ <j (x i -x j ), r 2 = -J2( x i-Xj) 2 , (8) 

i<j 
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and Pfc(x) is a homogeneous polynomial of degree k in the particle coordinates and 
satisfies the generalized Laplace equation i.e. 



N o 2 - N 
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Pk(x)=0 (9) 



As discussed in detail in [2], the polynomials Pk(x) are completely symmetrical under 
the exchange of any two coordinates. On inserting the ansatz (7) into the Schrodinger 
eq. (5) (with H given by eq. (1)) and using eq. (9) and following the procedure of ref 
[2] , we find that <f>(r) satisfies the equation 



-( " +S)0(r) = O (10) 
y JNr 



f(r) + {2k + 2b+l}-<P'(r) 

r 

where prime denotes differentiation with respect to the argument. The normalizable 
solutions of this equation are 

<P n ,k(r) = exp(-^[E~\r)L 2 n b+2k (2^[E~\r). (11) 

while the corresponding energies are as given by eq. (2). Here L%(r) is a Laguerre 
polynomial. Notice that in the expression (2) for the energy, n and k always come in 
the combination n+k (unlike in the Calogero case [2] where it comes in the combination 
2n+k). 

In the special case of N = 3, we can check our expressions for E n and VVi,fc with the 
exact expressions obtained by entirely different method (see eqs. (40) to (43) of [7] ). 
On comparing the two we find (note coupling constant in [7] is a/3« rather than a) that 
the two expressions agree provided k = 31 and Pk(x) oc r 3l C^ +1 ^ 2 (cos3(p) where Cf is a 
Gegenbauer polynomial. 

Let us now consider the positive energy spectrum of the Hamiltonian (1). It is of 
course purely continuous spectrum. Following the treatment given above and as in [2] 
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(Sec. 4) it is clear that the complete set of stationary eigenfunctions of the problem (in 
the center-of-mass frame) is 

V> pfc = Z a+1 / 2 <P p (r)P k (x), k = 0, 1, 2, p > (12) 

where p is connected to the energy eigenvalue by E = p 2 > (note that we have chosen 
h = 2m = 1) while <f> p (r) satisfies eq. (10). It is easily shown that for E > 0, the solution 
of eq. (10) is given by 

p (r) = e ipr F(k + b+ =, 2k + 2b + 1; -2ipr) (13) 

2 2p\/N 

One can now run through the arguments of [2] (Sec. 4) and show that if the stationary 
eigenfunction describing, in the center-of-mass frame, the scattering situation is charac- 
terized by the form 

N 

ip in ~ C(exp[i ^PiXi] (14) 
i=i 

with (note X{ > a^+i, i = 1, 2, N — 1) 



AT TV 



then -i/w is given by 



where 



i=i i=i 



AT 



Vw ~ Ce 2 ^-^exp[i ^pjv+i-iXi] (16) 



i=l 



,2i Vp _ 2 2pVN J 

" r(fc + 6+i + -3=) 



(17) 



2 2pVN> 

Thus we have the remarkable result that even in the presence of the potential 

, Q , the N-particle scattering problem reduces to a sequence of 2-body 

processes as characterized by eq. (4) but now one has an energy dependent phase shift. 



Note that all the results about scattering are also valid in case a is negative but now the 
spectrum is purely continuous and there are no bound states. 

Finally, let us discuss the "Sutherland variant" [3] of the Hamiltonian (1). Consider 

N d 2 a 

^"S^+IJf^F" vfe? (18> 

i.e. the N-body potential is now an external potential. Note that for the Calogero case, 
Sutherland [3] was able to obtain an exact expression for the ground state energy and 
eigenfunction ip and find a remarkable connection of ip 2 with the joint probability density 
function for the eigenvalues of matrices from a Gaussian ensembles in case (3 = 2A = 1, 
2 or 4. Using these connections he was the able to compute [3] the one particle density 
and the pair correlation function. 

Following Sutherland, let us consider the Schrodinger equation Hip = Eip with H as 
given by eq. (18). Further, let us write the wavefunction ifr as ip = 0<E> with 

(f) = n i<7 - | Xi - Xj | A , A = ^+a (19) 

On using this ansatz in the Schrodinger equation with H as given by eq.(18) we find that 
$ must satisfy 

d 2 $ _ 1 , d d 



It is easily verified that 



^ = exp(-y/\E\^^) (21) 



is a solution to eq. (20) with the energy 



e 4 



E [(iV-l)(l + AiV)]2 (22) 
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Clearly, for each ordering of particles, tfj is nodeless and hence it is the solution for the 
ground state. If we rewrite ip in terms of the variables 



where C is the normalization constant. A la original Sutherland case [3], where ip 2 was 
identical to the joint probability density function, it would indeed be remarkable if our 
ip 2 as given by eq. (24) for atleast (3 = 1, 2, 4 can be mapped on to some known solvable 
problem and using these results if one could obtain the one particle density and the pair 
correlation function for our case. 

This work raises several issues which need to be looked into. I hope to address some 
of these issues in the near future. 



Vi = 



Va 



(23) 



then one finds that 




(24) 
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